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_ Abstract 

We study the interaction of the a p and D*, D* with spins aligned using the Fixed Center Approximation 
to the Faddeev equations. We select a cluster of D*D* , which is found to be bound in / = 0 and can be 
associated to the X(3915), and let the p meson orbit around the D* and D*. In this case we find an / = 1 
state with mass around 4340 MeV and narrow width of about 50 MeV. We also investigate the case with a 
cluster of pD* and let the D* orbit around the system of the two states. The pD* cluster is also found to 
bind and leads to the £>2(2460) state. The addition of the extra D* produces further binding and we find, 
with admitted uncertainties, a state of I = 0 around 4000 MeV, and a less bound narrow state with 1 = 1 
around 4200 MeV. 

PACS numbers: 
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I. INTRODUCTION 


In hadron physics the existence of the three body resonances is drawing much attention for 
a long time. The proper analysis of the three hadron system can be tackled by combining the 
Faddeev equations [1] and chiral dynamics @-0]. However it is quite difficult to solve exactly the 
Faddeev equations. Recently, the combination of Faddeev equations and chiral dynamics was used 
to investigate for two meson-one baryon systems [s] and also for three mesons systems [(3, 3]. On 
the other hand, the fixed center approximation (FCA) to the Faddeev equations is technically 
very simple and a powerful method to explore three hadron systems. This method is especially 
suitable to study the system in which two of the three particles are bound forming a cluster and 
this cluster is not much altered by the collision of the third particle 0-0. This method has 


proved to be rather reliable for cases like A"-deuteron scattering very close to threshold [ill Il2j ]. 


In recent years the FCA to Faddeev equations has been successfully applied to the study of many 
three body interactions. For example in the work of the Ref. 13| the n — (A p)^ ri675) was 

analyzed by means of the FCA to the Faddeev equations in which the authors give a reasonable 
explanation for the A( 5 / 2 )+ (2000) puzzle. Likewise, the NKK system was investigated using the 


Faddeev equations under the FCA in Ref. and the results are in good agreement with the 


variational estimation [15] and also the full Faddeev calculation [la, [17(] . 

In the present paper we want to study the pD*D* system. The reason is that the vector-vector 
interaction is found to be very strong, particularly in the J = 2 channel [181 . [la]. Due to this, it 
was possible to see that multi-p states with the spins parallel were bound, although the width was 
increasing with the number of p’s {2d]. The masses and widths obtained were in good agreement 
with experimental data. The same occurred with K* multi-p states in |2l( and with D* multi-p 
states in |22|. In the case of K* multi-p states one finds good agreement with experiment, but in 
the case of D* multi-p states only predictions were made and experimental counterparts have not 
yet been observed. 

The pD*D* system is new and has not been studied so far. Yet, studies done for D*D* 
interaction in [23| and D*p interaction in [2J] have already set the grounds to tackle this interesting 
system with hidden charm, and we wish to study it here. 


II. FORMALISM 


In this section we describe the calculation of the three body interaction of the pD*D* system 
in s-wave and all the spins aligned. In the work of [241 ]. the pD* interaction was studied using the 
hidden gauge formalism [25H27|. In 2J] the authors found strong attraction in / = 1/2, J = 2 
which corresponds to the tensor state A>2(2460) with I(J P ) = ^(2 + ). Similarly, it was also found 
in Ref. Q that the resonance Y(3915) ( I G (J PC ) = 0 + (2+ + )) could be understood as a molecule 
made of D* and D* mesons with strong attraction in the I = 0, J = 2 sector. 

We use the FCA to the Faddeev equations to study the pD*D* system. This method is partic¬ 
ularly well suited for the system in which a pair of particles clusters together and the cluster is not 
much modified by the third particle, like in the case of D2(2460) as a pD* cluster and 2f(3915) as 
a D*D* cluster. 

For the technical details we proceed similarly to Ref. [ld| to express the FCA to the Faddeev 
equations. Hereinafter we are going to interpret the interaction of a particle a 3 with a cluster made 
of two particles, ai and 02 - In this work the particle 03 will represent the p(D*) which scatters 
from the cluster, X(3915) (A>|(2460)) and ai and <22 are the D*(p ) and D*(D*) which build up 
the cluster. 

The diagramatic representation of the FCA to the Faddeev equations is shown in Fig. [T] The 
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FIG. 1: Diagrammatic representation of the fixed center approximation to Faddeev equations. 


particle 03 rescatters repeatedly with the components of the cluster. In Fig. [T| the thick squared 
dots represent the unitary scattering amplitudes with coupled channel for the interaction of particle 
<23 with particle ai (ti) and 02 (£ 2)1 respectively, which will be discussed later. In order to write the 
equations for the total three body scattering amplitude, we define two partition functions T±, T 2 
which sum all diagrams of the series of Fig. [I] which begin with the interaction of particle <23 with 
particle ai of the cluster (Ti), or with particle (22 {T 2 ). Then the FCA equations can be written as 
a system of coupled equations: 

T\ = fi+tiGoT^, 

T 2 = t 2 + t 2 G 0 T u 

T = Ti + T 2 (1) 

where Go is the Green function for the propagator of particle <23 between the particles ai and a 2 
which is discussed later on. 


The scattering amplitude (pD*D*\t\pD*D*) for the single scattering contribution is obtained in 
terms of the two-body amplitudes t \, t 2 derived in Refs. 23, 24]. First we explicitly determine the 
case of p(D*D*) in which = 0 and the total isospin of the three body system I p ( D *£,*) = 1. 

Using the nomenclature | p,I,I z )(g> | DD,I,I Z ) we obtain 
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with the nomenclature | I z ) for the p meson and | I Zl ,I Z2 ) for the D*D* system. 

Second we write the case of D*(pD*) where I p d* =1/2 and the total isospin of the three body 
system If)^ pD *) = 0 or I^( p n*) = 1 - 
For the total isospin 1 = 1 case 


(D*(pD*)\t\D*(pD*)) = (D*(pD *) I=1/2 | (tVp + ^W) I D*(pD*) I=1 / 2 ) 
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with the nomenclature | D *, I z 
similarly derive 


pD , I Z1 , I z 2 ) ■ In the case of the total isospin 7 = 0 case we 


(D*(pD*)\t\D*(pD*)) = (ty/ 2 ) + (^). 


( 5 ) 


Since we use the normalization of Mandl and Shaw [28j], which has different weight factors for 


the particle fields, we need to consider how these factors are adapted to the present problem. It is 
easy to do this compairing the single scattering, double scattering and full scattering amplitudes. 
In this case, following the field normalization of Ref. [28l ]. we can obtain the S matrix for the single 
scattering diagram (Fig. [Q (a) and (e)), 


5*7 = -iti^(2nYS 4 (k + k R - k' - k' R ) 

1111 


( 6 ) 
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~ it2 ^2 ( 27 r ) 4 <5 4 0 + k R - k' - k ' R ) 

1111 


where the momentum k(k'), the on-shell energy refer to the initial (final) particles, respec¬ 

tively, and V is the volume of the box where the states are normalized to unity. In Eqs. © and 
m, t \, £2 correspond to the first and second terms of the right hand side of Eqs. 0 and (f5|). 
Likewise we have the 5-matrix for the double scattering diagram as (Fig. [T] ( 6 ) or (/)) 


5 (2) = -i(2ir) 4 ^5 4 (k + k R -k' -k' R ) 

111111 

/ d 3 q _, . . 1 

x / —^2- 2 , • 2, 

J (Z7TJ- 3 —q^ — m^+ie 


( 8 ) 


where F R (q ) is the form factor of the cluster which represents essentially the Fourier transform of 
its wave function. The derivation of the form factors can proceed similarly to Refs. 0,0 , where 
one can also read further discussions and interpretation. The form factor for s-wave functions is 
given by 


F R (q) = Tr / P<kmax 

■N J\p—$\<kma 

1 


d 3 p 


1 


1 


1 


2u ai (P ) 2 U>a 2 (p) M R - UJ ai (p ) - UJa 2 (p ) 

1 \ 1 


,2 u] ai (p-q)J \2 u a2 (p-q)J M R - u 01 (p-q) - uj a2 (p~ q) ’ 
with the normalization J\f 


(9) 


M = 


Jp<k„ 


d 3 p 


1 


1 


1 


2u ai {p) 2UJa 2 (p) M R -UJ ai (p )-UJ a2 (p). 


( 10 ) 


where u> ai and uj a2 ar e the energies of the particles a\, 02 , and k max is a cutoff that regularizes the 
integral of Eqs. @ and (1101) . This cutoff is the same one needed in the regularization of the loop 
function of the two particle propagators in the study of the interaction of the two particles of the 
cluster 01 . In this work we take the cutoff k max = 1200 MeV, the same one used to generate the 
Z ?2 (246oj [24 ] . 

Similarly, the full three body 5-matrix for scattering of particle a 3 with the cluster is given by 


5 = -iT^(2ir) 4 5\k + k R -k' -k' R ) 


1 


1 


1 


1 




(ii) 


Compairing this equation with Eqs. 
amplitudes: 


and 0 , we introduce convenient factors in the elementary 
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with m ai , rn a2 and M R the masses of the particle a\, a 2 and the cluster respectively, where we 
have taken the approximations, suitable for bound states, 


/ 


2 cjg 


1 ( 2 ) 


V 


2 m a 


1 ( 2 ) 


Finally solving the set of equations for the FCA to the Faddeev equations, Eqs. (pQ) , we obtain 


T = T 1 +T 2 


+ t 2 + 2t\t 2 Go 

1 — t\t 2 G'y 


(13) 


Note that the argument of the total amplitude T is regarded as a function of the total invariant 
mass of the three body system whereas the argument of 2 ) is the invariant masses of the two 
body systems. In order to obtain the arguments s^) of the two body amplitude we share the 
binding energy among the three particles, proportionally to their masses. Therefore the energy of 
the particles ai, a 2 and 03 become 


F -m ^ 

^- m a^M R + m a3 ) 


(14) 




m ai M r 


E = 

ai (M R + m a3 ) (m ai + m a2 ) 


E = \fs _ m a2 M r 

“ 2 (M R + m a3 ) (m ai + m a2 ) 


Hence the total energy of the two body system is evaluated as follows 
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with H a2 ( ai ) the binding energy of the particle a 2 (ai). 

As we stated before, the Go function is the propagator of the particle 03 inside the cluster as 
follows 


Go 


1 

2 M r 



(19) 


where M R is the mass of the cluster, and m a3 the mass of the particle 03 . Here the energy of the 
propagator q° is determined at the three body rest frame as 


q 


0 


s + m 


a 3 


-M 2 r 


2 Vi 


( 20 ) 
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FIG. 2: Real (solid line) and imaginary (dashed line) parts of Go function in p(D*D*) 


with y/s the rest energy of the three body system. As an example we depict in Fig. [2] the real and 
imaginary parts of the Go function for the p(D*D*) system. The Go function has a similar shape 
for the D*(pD*) system but with a different threshold. 


As we mentioned previously, in the evaluation of the pD*D* three body interaction, the use 
of the pD*, pD* and D*D* unitarized amplitudes has crucial importance. These amplitudes were 


studied by the coupled channel Bethe-Salpeter equations in Refs. [23|, |24(] and we use them here. 


In order to reproduce the pD*(pD*) system from the work of [241 ]. the coupled channels used are 


pD* and uD* for I = \ and pD* for / = | case. In the case of the D*D* system [19j, [23[ there are 
10 coupled channels, D*D *, K*K*, pp, ww, (jxp, J/TJ/T, uJ/\V, ucj) and D*D* for 1 = 0 

and six coupled channels D*D *, K*K*, pp, pu, pJ /T and pcf) for 1 = 0. 

Following the ideas of the coupled channels chiral unitary approach, the FF-two body scattering 
amplitude can be obtained using the Bethe-Salpeter equations in its on-shell factorized form as 
below 


t = (1 — FG) -1 F 


( 21 ) 


where the F is a matrix of the interaction potentials between the channels, which is calculated 
from the hidden gauge Lagrangian 251-27]. The potential F is a 10 x 10 matrix in / = 0 and 6x6 
matrix in I = 1 with the amplitudes obtained from the coupled channels for D*D* case in J = 2 
0, 0]. In addition, in the case of the pD* the potential F is a 2 x 2 matrix given by [24| . 

In Eq. (I21|) G is a diagonal matrix of the loop function of two mesons in the i channel 


Gi(P) = i 


d A q 1 1 

(27t) 4 q 2 — mi + ie (P — q) 2 — m 2 + ie 


( 22 ) 
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where P is the four dimensional momentum of the two vector mesons determined at the rest 
frame, P = (yfs, 0), and mi and m 2 are the masses of the vector mesons in the i channel. In order 
to remove the ultraviolet divergence of the loop function, we use the dimensional regularization 
scheme and we get 
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(23) 


where qi is the three momentum of the two vector mesons determined at the center of mass frame 
evaluated as follows 


Qi = 


a/(s — (to-i + m 2 ) 2 ) (s — (mi - m 2 ) 2 ) 

Vi ' 


(24) 


In Eq. (1231) . g is a regularization scale and a, is the subtraction constant. We take g = 1500 
MeV and a = —1.74 for the pD* to reproduce TAJ (2460) resonance 24] and p = 1000 MeV and 
a = —2.07 in the channel including heavy mesons, a = —1.65 in the channel including light mesons 
for D*D* to obtain the X(3915) state [23|]. In Eq. @, the cutoff regularizes the integral of the 
form factor. In the present work, we choose the cutoff k max = 1200 MeV both for the L>2(2460) 
and X(3915), which produces similar results as the using the chosen subtraction constants in the 
dimensional regularization. 


III. RESULTS 


In this section we present the results obtained for the scattering amplitude of the pD*D* system 
in spin-3. The two-body pD*, pD* and D*D* systems were investigated by the coupled channel 


Bethe-Salpeter equations in Refs. [231 . [24j]. As we stated before, the resonance D%( 2460) was 


generated as a pD* quasibound state or molecule in the isospin 1/2 and spin-2. It was also found 
that the resonance X(3915) is dynamically generated in / = 0 and spin-2 from D*D* scattering. 
Therefore, there are two possible cases of three-body scattering of the pD*D* system. One is the 
IAj(2460) — D* and the other one is the V(3915) — p. 

In Fig. [3] we illustrate the modulus squared |T| 2 for the Xp —»• Xp scattering as a function of 
the total energy of the pD*D* system for the case of / = 1 and J = 3. The results show a clear 
peak at y/s = 4338 MeV about 360 MeV below the threshold of the A(3915) — p system. The 
width of the peak is about 50 MeV. 

This looks like a strong binding, but we must keep in mind that the vector-vector interaction in 


J = 2 is indeed very strong [18|, |l9|, [23|, 24(. This is why we are studying these superbound states 


with spins aligned where the spin of any pair is always J = 2. The obvious thing is that the D*D* 
state is already bound and since the pD* also binds to give the I?|(2460), the system pD*D* with 
this configuration will necessarily be bound. This would be the case even if the p interacted only 
with one D*. In this case we would have a binding of mo* + wip — mp* = 320 MeV. The binding 
that we get is 360 MeV, which means that we have gained extra 40 MeV binding by the interaction 
of the p with the D*. This indicates that there is extra binding from the three body molecular 
structure. This feature is reminiscent of what happens in Quantum mechanics for the problem of a 
particle in a well of two attractive (5-functions [30]. For the symmetric solutions, if the two <5’s are 
separated, the binding energy of the two 5 potential is the binding of a particle in one <5 well. As 
















FIG. 3: Modulus squared of the p(D*D*) scattering amplitude with total isospin 7 = 1. 


the two 5 potentials get closer, the particle starts orbiting the two potential wells and the binding 
energy grows. This is what we observe here, indicating an extra binding from the orbiting of the 
p around the D* and 77*. 


In Figs. [4] and [5] we show the results of |T| 2 for the case of 77|77* —> D* 2 D* in total isospin 
7 = 0 and 7 = 1, respectively. We find a peak around 4000 MeV which is about 470 MeV below the 
7?2 (2460) and 77* threshold for the isospin 7 = 0 case. The width of this state is quite large about 
250 MeV. If we conduct the same exercise as before, the 77* with a p would be bound by 320 MeV, 
and the 77* with the 77* by about 63 MeV. We would think that the 77* is orbiting the p where 
it is more bound and get some extra binding from orbiting the 77*. It is not clear why one passes 
to 470 MeV binding. It is also unclear why the width is much larger. Probably we have to accept 
that in this case, since the 77* is heavier than the p, there are limitations to the applications of the 
FCA, and we should accept this result as an indication that we could now have a state more bound 
than in the former case, which fulfills all the conditions for a reliable application of the FCA and 
hence is more reliable, but we cannot be certain about the mass and the width. 

For the isospin 7 = 1 case, we see a clear peak around 4195 MeV, and the width is around 
60 MeV. The position of the peak is about 270 MeV below the 77|(2460) and 77* threshold. 
These results are more intuitive than before. We can apply the same argumentation as before, 
but now according to Eqs. (JJ|) and <[5|), we can see that the weight of the amplitude in 

(77*(p77*)|t|77*(p77*)} for I = 0 is unity while for 7 = 1 it is 1/9, and this is the amplitude that 
contains the attractions that binds the 772(2460). With the caveat about the arguments used before 
for the case of 7 = 0, it looks clear that the binding should be smaller than for the case of 7 = 0 
and the width is also similar to that of the p(77*77*) molecule. 
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FIG. 4: Modulus squared of the D*(pD*) scattering amplitude with total isospin 1 = 0. 



FIG. 5: Modulus squared of the D*(pD*) scattering amplitude with total isospin 7 = 1. 
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IV. CONCLUSIONS 


We studied the interaction of the p and D*, D* with spins aligned using the Fixed Center 
Approximation to the Faddeev equations. This guarantees that we have J = 2 for any of the pairs, 
where the vector-vector interaction is strongest, and leads to three body states with J = 3. We 
first select a cluster of D*D*, which is found to be bound in / = 0 and can be associated to the 
X(3915), and then let the p meson orbit around the D* and D*. In this case the FCA produces an 
amplitude for p-D*D* scattering which has a clear and narrow peak around 4340 MeV. The case of 
a D* orbiting around a cluster of pD* is more uncertain because the mass of the external particle 
is heavier than the one of the p in the cluster, and the FCA is less reliable. In this case the cluster 
makes the 2460) state, and we point at some qualitative results, with an I = 0 state around 
4000 MeV and an I = 1 state around 4200 MeV. The results obtained for the 1 = 1 state with the 
p orbiting around the X(3915) should be realistic since the p is lighter than the constituents of the 
cluster. In the other case our results should be taken as indicative, but strong arguments are given 
that these states should be strongly bound. 

The results obtained here should serve to stimulate calculations with more accurate three body 
tools, as those of which could make predictions on these interesting J = 3 states. One 

should recall at this point that states with increasing spin number already exist in the light sector 
[2o| and in the strange sector [2l|. What we have done here it to extend this to the hidden charm 
sector. Parallelly, it would also be interesting to investigate states of large spin in the region of 
mass investigated here. The results obtained in this work provide sufficient support for a devoted 
search of such states. 
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